Boundary conditions for spin diffusion 
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We develop a general scheme of deriving boundary conditions for spin-charge coupled transport 
in disordered systems with spin-orbit interactions. To illustrate the application of the method, we 
explicitly derive boundary conditions for spin diffusion in the Rashba model. Due to the surface 
spin precession, the boundary conditions are non-trivial and contain terms, which couple different 
components of the spin density. We argue that boundary conditions and the corresponding electric- 
field-induced spin accumulation generally depend on the nature of the boundary and therefore the 
spin Hall effect in a spin-orbit coupled system can be viewed as a non-universal edge phenomenon. 
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I. INTRODUCTION 



There has been a lot of recent interest in a phenomenon 
dubbed the spin Hall effeclii* 2 * 3 ". Experimentally this ef- 
fect manifests itself as an equilibrium spin accumulation 
near the edges, when an electric field is applied paral- 
lel to a Hall bar of a spin-orbit coupled electron sys- 
tem. This kind of electric-field-induced spin polarization 
has been observed in recent experiments^. This phe- 
nomenon has a potential for being technologically im- 
portant, as it allows one to electrically manipulate the 
spin degree of freedom in non-magnetic systems. There 
are still many controversies regarding the physical ori- 
gin of the spin accumulation observed in real experiment 
(intrinsic 1 " vs. extrinsic 2 '^). But even within a given the- 
oretical model (e.g., a two-dimensional electron gas with 
Rashba or Dresselhaus spin-orbit interactions or Lut- 
tinger model), there is still no consensus on the correct 
description of the spin Hall effect 7 ". 

A common theoretical approach to the intrinsic spin 
Hall effect is as follows: One introduces the notion of a 
spin current, which is an operator intuitively related to 
spin transport. Using the Kubo formula, one then calcu- 
lates the spin-current autocorrelator and the correspond- 
ing linear response quantity, called the "spin Hall conduc- 
tivity'— In clean systems, it appears to be related to an 
elegant topological Berry's phase structured. In a dis- 
ordered system, the notion of the Fermi surface Berry's 
phase is unclear, but the "spin Hall conductivity" is well- 
defined. It is believed that this quantity is connected to 
spin accumulation in the spin-Hall experiment. While 
this theoretical approach is mathematically well-defined 
and may lead to technically challenging problems^, its 
relation to experiment remains unclear. The problem is 
that in spin-orbit coupled systems, the spin is not con- 
served and thus the definition of the spin current is am- 
biguousiSiii. A non-equilibrium spin density in a region 
relaxes not only due to the flux through the boundaries 
of the region, but also due to spin precession. Therefore, 
the spin current is not easily measurable. 



To describe realistic experiment, it is preferable to use 
another route without involving the notion of a spin cur- 
rent. In disordered systems, this can be done with the 
help of a kinetic equation^ 2 - or a spin-charge coupled dif- 
fusion equation 13-14 . The diffusion equation provides a 
complete and physically meaningful description of spin 
and charge diffusive transport in terms of position and 
time dependent spin and charge densities. However, the 
explicit solution of the diffusion equation (which is gen- 
erally a complicated set of partial differential equations) 
strongly depends on boundary conditions 15 . While for 
usual charge transport, boundary conditions are often ob- 
vious and follow from conservation laws, this is not the 
case in spin-orbit coupled systems. The problem is that 
the spin is not conserved and thus there is spin preces- 
sion at the edges. This spin dynamics is very sensitive 
to the actual boundary (depending on the physical sit- 
uation, the boundary could be modelled as a hard wall, 
WKB potential, rough surface, etc.). Since the bulk dif- 
fusion equation describes spin/charge dynamics only at 
the length-scales much larger than the mean- free path I, 
it can not capture the behavior of the spin in the imme- 
diate vicinity of the sample edges. Due to the aforemen- 
tioned boundary spin precession, the "trivial boundary 
conditions" (i.e., boundary conditions, which one would 
have, if the spin had been conserved) acquire corrections 
proportional to the spin-orbit coupling. These correc- 
tions are very important, as they select a unique solution 
of the diffusion equation and determine the observable 
spin accumulation. 

The derivation of boundary conditions is technically 
a difficult task. Generally, one needs to determine the 
behavior of the system in the ballistic region near the 
boundary and match it with the diffusive behavior in the 
bulk. It requires the knowledge of the exact S'-matrix 
or the boundary Green's function. We should point out 
that there exists an extensive literature on the related 
issues in the context of the diffusion of light in a ran- 
dom medium 1 ^ and neutron diffusion 17 . In these prob- 
lems, the general form of boundary conditions is obvious 
due to conservation laws and only numerical coefficients 
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are unknown. There exists a general method of deriving 
these numerical coefficients. But there are just a very 
few models in which exact analytical results are avail- 
able (a notable example is the Milne problem 16 , which 
describes a boundary separating diffusive and ballistic 
regions). We note that in the context of spin diffusion 
in spin-orbit coupled systems, even the qualitative form 
of boundary conditions is not known as there are no ob- 
vious conservation laws, which would provide guidance. 
Even though in the bulk, one has conserved quantities 
labeled by the band index (e.g., chirality in the Rashba 
model), the boundary generally does not respect these 
conservation laws and mixes up different bands. 

In this paper we formulate a general method of deriving 
boundary conditions in spin-orbit coupled systems. Us- 
ing this proposed method, we derive boundary conditions 
for the Rashba modeli^ in the leading order with respect 
to the spin-orbit coupling, which is assumed small. We 
show that there are corrections to the "trivial bound- 
ary conditions:" In leading order, we find that at the 
boundary, the spin component perpendicular to the plane 
of the two-dimensional gas gets coupled to the in-plane 
spin component perpendicular to the boundary. We ar- 
gue that spin accumulation in the spin Hall effect setup 
is determined by the combination of the bulk spin-charge 
coupling and boundary effects. 



II. GENERAL METHOD OF DERIVING 
BOUNDARY CONDITIONS 



In this section we formulate a general method of de- 
riving boundary conditions applicable to any spin-orbit 
coupled system with a boundary. Let us consider a dis- 
ordered electron system occupying a region A. One can 
derive the following integral equation for the diffuson 19 , 
which determines the dynamics of the charge density (po) 
and the spin densities (p a , a — x,y, z): 

X> a b(w;ri,r 2 ) / V cb (uj;r ,r 2 )d a r 



u r'eA 

-S ab S{r 1 - r 2 ) 



(1) 



where v is the density of states per spin at the Fermi 
surface, r is the scattering time, the Latin indices label 
charge (0) and spin (x, y, and z) degrees of freedom, 
and T> a b(co; r%, r 2 ) is the diffuson, which is the Green's 
function of the spin-charge coupled diffusion equation. 
The kernel in Eq. Q reads 

n afc ( W ;r,r') = ^Tr {a a d R (e 1 ;r,r')a b d A (e 2 ;r' ,r)} , 

(2) 

where lo = E\ — e 2 , &o is the unity matrix, a a are the 

Pauli matrices, and G R / A (e; r, r') are disorder-averaged 
retarded/advanced electron Green's functions, which in- 
clude the effect of the boundary (here and below we de- 
note the boundary Green's functions as G and the bulk 



Green's functions as G). They depend exponentially on 
distances via factors of the following types, cx e~ r ^ 2l \ 
where I is the mean-free path. Deep in the bulk, the 
edge effects are exponentially small and the bulk Green's 
functions can be used, G R I A (e, |r — r'|). Performing a 
gradient expansion of the bulk polarizability (i.e., tak- 
ing the limit I — * 0) , one obtains the spin-charge coupled 
diffusion equation. Generally, it has the following form: 



d tPa = (D a d 2 -T- 1 )p a + J2ri Y d. 



'iPb, 



(3) 
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where p a (r,t) are the charge/spin densities, D a are the 
diffusion coefficients, r a are the relaxation times, and T^, 
are the spin-spin and spin-charge couplings. 

Note that the integral equation Q does not require 
boundary conditions (it already contains this information 
through the boundary Green's function), but the differ- 
ential diffusion equation docs. To derive these bound- 
ary conditions one needs to perform the gradient expan- 
sion near the boundary dA. Let us consider the problem 
in half-space A — {x \ x > 0}. To find the diffusion equa- 
tion near the boundary, we take the limit x = r • n — > 0. 
This limit should be understood in the following sense: 
Pp 1 « i < I, i.e., we consider a point, which is near the 
boundary as compared to the diffusion length-scale I, but 
still far from it as compared to the ballistic length scale 
described by the Fermi wave-length. Let us first consider 
the case of zero spin-orbit interactions. If the boundary 
is a hard wall potential, then the boundary Green's func- 
tion can be obtained using the method of mirror images 



6 R/A (e; r, r') = \G^ (e; Ar) - Gq'" (e; Ar,)j <7„, 

(4) 

where Ar = r — r', Ar* = r* — r', and is the mir- 
ror image with respect to the sample boundary. For the 
sake of concreteness, let us assume that the dimension- 
ality d = 2 (all qualitative arguments are independent of 
dimensionality). In two dimensions, the Green's function 
has the following form 



yR/A 



R/A 



G R ' A (e;r) = -, 



' ztim 2 



exp 



±i \Pf H 



(5) 

We note here that the frequency dependence of the dif- 
fusion kernel is irrelevant for the boundary condition 
problem and will be omitted from now on. When the 
method of mirror images applies, the polarizability kernel 
(J2J) has the following general structure II^ ) = G R G A + 
G R G A —G R G A —G R G A [where G* implies the "reflected" 
Green's function, i.e, the second term in Eq. The 
key observation is that due to the fast-oscillating fac- 
tor GG* oc exp [±ipF (Ar — Ar*)] and our choice of the 
point r (x pp 1 ), the cross terms GG* in the expression 
for IT ) become negligible in the integral over r' and can 
be omitted. The remaining integral can be easily calcu- 
lated and we find (in the "trivial" case of no spin-orbit 
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coupling) : 

li m V / ° 6 ' Pb{T')d d v' = p a (0) + c d l d xPa (0), 

b x>>0 

(6) 

where a is a dimensionless coefficient, which in two di- 
mensions is equal to c-i — 2/ir. From Eqs. and JBJl, 
we recover the familiar boundary condition n -dp(0) = 0, 
which is simply the consequence of the charge and spin 
conservation. 

In a system with spin-orbit interaction, the method of 
images does not apply because the "reflected" Green's 
function does not satisfy the corresponding Schrodinger 
equation. This is a major complication, which never oc- 
curs in the problems usually studied in the literature in 
this contex1ji£ii£. In addition, the matrix structure of 
the boundary diffusion kernel becomes essential. There 
are generally two ways to obtain the boundary Green's 
function (we have explicitly verified that both approaches 
lead to identical perturbative results in the leading or- 
der—): (i) The first is to solve the Schrodinger equation 
near the boundary and express the Green's function using 
the corresponding wave-functions as a basis; (ii) The sec- 
ond is to find the Green's functions perturbatively, using 
the mirror-image result Q as the leading approximation: 

<5G(r,r')=]T I ■ I d d r 1 ...d d r k G (r,r 1 )V so (r 1 ) 

k= \: k >0 x t >0 

xG (ri,r 2 ) x ... x Vs O (r fc )G (r fe , r'), (7) 

where Vso is the spin-orbit interaction operator as it ap- 
pears in the Hamiltonian. We note that disorder correc- 
tions to the spin-orbit vertex can be shown to be small 
as long as EpT ^> 1. The resulting Green's function 
Q may be quite complicated even in the simplest cases. 
However, the diffusion kernel always contains products 
of the following types: 

n afc (r, r') oc [TGq (Ar) - iH?£(Ar„)] 

x [T'G*(Ar)-R'G£(Ar*)], (8) 

where T's and i?'s are some functions, which do not oscil- 
late on a Fermi wavelength scale. Again, in the integral 
of Eq. QJ, the cross-terms G^ A (Ar)G^ R (Ar*) become 
negligible due to the fast-oscillations and can be omitted, 
but the other terms survive. One can define the following 
matrix, which determines boundary conditions 

„ ,. f SU a b(r,r') ld , . . 

B ab = hm / — dV, 9 

r^O J 2-KVT 

x'>0 

where the correction to polarizability is determined by 
the Green's function Q and Eq. and the limit is 
understood in the sense pp 1 <C |r-n| <C I. If the boundary 
is an impenetrable wall, general boundary conditions take 



the form 

n-dp a (0) = -^-J2 B ^Pb{0), (10) 
c d i b 

where a,b = 0,x,y, and z, n is a unit vector normal 
to the boundary, is a number (c 2 = 2/7r), and B a \, 
is a dimensionless matrix defined above. If an external 
field is present one may need to consider another term in 
the boundary conditions, which has the meaning of the 
boundary spin-charge coupling 21 

C =lmi / <m ah (r,r')r'dV/(27rz/T). (11) 

r^O J 

x'>0 

In this case the boundary condition becomes: 

c d ln ■ dp a {0) =-J2 B ab Pb(0) - C Q Vp (0) (12) 

b 

We note that if the boundary is not impenetrable, 
the general method does not change but the results do 
even in the lowest order. For instance, if the bound- 
ary is characterized by an isotropic reflection coefficient 
R =/= I, the zeroth order boundary condition becomes 
(1 - R 2 )p a (0) = c d l(l + R 2 )n ■ dp a {0). If the reflection 
coefficient is not too close to unity, the right-hand side of 
the latter equation can be considered small and one can 
use the following leading order boundary conditions 

Pa(0) = 0. (13) 

This boundary condition implies that an excess charge or 
spin density can not exist at the boundary if the bound- 
ary is penetrable. The excess density leaks through the 
edge (see also Ref. |22j). 

We note that strictly speaking the matching point at 
the boundary [which we set to zero in Eqs. jjEt and fT3)l ] 
can not be obtained exactly within the method described 
here, but all dimensional parameters and numerical co- 
efficients are robust. To determine the correct matching 
point on a lengthscale of order I, one actually needs to 
solve exactly the integral equation, which is possible only 
in very special cases [e.g., the Milne problem (i.e., K—0 
and no spin-orbit coupling)]. 



III. APPLICATION TO THE RASHBA MODEL 

To illustrate the application of our method, let us con- 
sider the Rashba modeU^ in the leading order in spin- 
orbit coupling. In this case, the spin-orbit interaction 
operator has the form Vso = —ict^ie a i3 z a a dfs, where 
is the Rashba spin-orbit coupling and a, (3 = x, UjZ. The 
differential diffusion equation was derived in Ref. |l3l | and 
reads 

d t po = Dd 2 p - T cs e zat3 d a pp; (14) 
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d t p a = {Dd 



T a 1 ) Pa+^ss^ap-i^zXpdxp^ — T sc e z f} a d p Pq. 

(15) 

The spin-spin coupling in the Rashba model is r ss = 
Aa-^E-pT and the spin relaxation times are 2t z = t x = 
T y = r s = 2r/(2aRPFT) 2 . These parameters are univer- 
sal in the sense that they do not depend on the details 
of the electronic spectrum and can be calculated in the 
^-approximation^. However, the spin-charge coupling 
is not universal, since it strongly depends on the elec- 
tronic spectrum and the high-energy cut-off (if the spec- 
trum is quadratic, r sc = 2a R p|r 2 ). We should point out 
however that this result holds only in the semiclassical 
diffusion approximation, i.e., only in the leading order 
with respect to the inverse conductance 1/(Eft), which 
is an independent parameter of the model. There may 
exist contributions to the spin-charge coupling of order 
r sc ~ T ss /(Eft). They originate from diagrams with 
crossed impurity lines, which are omitted in the diffusion 
approximation (we note that there is no direct relation 
between r sc and the spin Hall conductivity, which van- 
ishes to all orders^) . These quantum corrections can be 

neglected only if {Eft) ' -C apFT <C 1. This is a very 
strong constraint, which is not always satisfied in exper- 
imentally studied systems. In what follows, we will treat 
the spin-charge coupling terms in the diffusion equation 
(|14l I15f) as non-universal phenomcnological parameters 
and assume that the above-mentioned constraint is sat- 
isfied. 

The diffusion equation iflil H3)l must be supplemented 
by boundary conditions. We consider the problem in 
leading approximation with respect to the spin-orbit cou- 
pling parameter ckr and the inverse conductance. First, 
we use Eqs. (01 and (0) to derive a correction to the po- 
larizability kernel 

<&n$(r, r') = -2a R e Q/W Re [<^(r, r')g 5 {r, r')] (16) 

and 

OTg(r,r') = 5ng(r,r') = 2a R ep S lm [g$(t, r')g s (r, r') 

where the boundary Green's functions are defined by 
Eqs. and © and we introduced the following function 



d l T 1 Gif(T,T 1 )d la Gg(r 1 ,T>). 



This function contains four terms (we can symbolically 
denote them as gi a — J Gd a G, g2a = — J G*d a G, 
g-3a = — J Gd a G„, and g^ a = J G*d a G*). Let us an- 
alyze the first one, the other terms can be calculated in a 
similar manner: To calculate the corresponding integral, 
we use a new coordinate system (x, y) in which the points 
r and r' have the coordinates (0, — Ar/2) and (0, Ar/2) 
respectively (Ar = r — r'). Next, we introduce the ellip- 
tic coordinates (77, </>), with the points r and r' in the foci 



of the ellipses: |ri — r| + |r' — ri | = i]Ar and <j> being the 
angle between the vector ri and the x axis. Using these 
elliptic coordinates and Eq. 10, we write the integral gi a 
in the following form 



5ia(r,r ) = -- d a 

AirpF 



Ar J dr] 

1 X ± >0 

x e 



/ r/ 2 — sin 2 1 
if - 1 



(18) 



We note that the limits of integration over <p are generally 
non-trivial due to the constraint x\ > 0. It is possible 
to evaluate analytically the correction to the diffusion 
kernel [see Eqs. ©, {SJ, and l|17|)] and the boundary 
matrix ©. Indeed, from Eqs. @, JEJ, and ifTTjl. we 
find that due to the fast-oscillating exponents, only a 
few terms survive the integration over r'. These remain- 
ing terms contain combinations of the following types 
exp [ippAr(?7 — 1)], which constrain the parameter 77 in 
(|18fl to be close to unity. We find in the leading order 
the following results for the boundary matrix |[SJ| 



4a R mZ 
-Dq/3 = e 



uj3y 



and Bn h = B h n = Sh„0 



E f t 



(19) 

where the last equation implies that the coupling between 
the spin and charge densities is small with respect to the 
inverse conductance and vanishes in the framework of the 
diffusion approximation and the accuracy of the method. 

From Eqs. IjlUfl and (|19|) . we find the following bound- 
ary conditions in the leading approximation with respect 
to the spin-orbit coupling and the inverse conductance: 



d xPo {Q) = d xPy (0) = 0, 



(20) 



Wo) = -§jp*(°)> and d *pM = ( 21 ) 

where r ss is defined after Eq. (I15|) . These boundary con- 
ditions can be generalized to describe any form of a hard- 
wall boundary. To find such a general form of bound- 
ary conditions, one should do the following replacements: 
d x — > n ■ Q and p x — > p n , where n is the unit vector nor- 
mal to the segment of the boundary and p n is the corre- 
sponding component of the spin density. However, if the 
boundary is a smooth potential with the length-scale in 
the n-direction being a, then boundary conditions would 
be different and strongly dependent on the ratios of a/ 1 
and app^S. Also, the boundary potential would lead to 
an additional boundary spin-orbit coupling, which would 
strongly affect the boundary conditions. 

We note here that in the next-to-leading order with 
respect to the spin-orbit interaction, two new effects ap- 
pear: Boundary spin-charge coupling and boundary spin 
relaxation (the former is parametrically larger than the 
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latter due to an additional smallness of the spin density) . 
To this order boundary conditions take the form: 



and for the hard-wall boundary 



n-d Pz (0) 



2D 



P„(0)+7P,(0)-C z Vp (0), 



where 7 and C z are constants oc a^. We should note 
that the boundary spin-charge coupling term can be 
shown to be strongly dependent on the electronic spec- 
trum and the physics far from the Fermi line and as 
such is non-universal and can not be calculated in the 
^-approximation used in this paper. 



IV. SOLUTION OF THE DIFFUSION 
EQUATION 

We now solve the diffusion equation for a two- 
dimensional gas of Rashba electrons occupying a half- 
space x > and in the presence of an electric field par- 
allel to the edge. We consider two types of edges: a 
partially or totally transparent interface and an impen- 
etrable boundary. It is convenient to rewrite the dif- 
fusion equation (|14l I15f) in terms of the rescaled vari- 
ables: p x (x) = [ p x (x) - poo]/ poo and p z (x) = p z (x)/p QOl 
where x, = xj 1 \j ' Dt s is the dimensionless distance and 
Poo = 2T sc vT s eE is the uniform spin density due to the 
spin-charge coupling. In these variables, the diffusion 
equation takes the following form 



Px 



2p' z = 0; 



f z - 2p z - 2p' x = 0. 



(22) 



(23) 



Boundary conditions for a transparent boundary are 
p z (0) = and p x (0) = — 1 and for a hard- wall bound- 
ary ^(0) = -p z {0) and p'M = p x (0) + C, where C = 
1 — C Z ydyPo/ poo is a non-universal spectrum-dependent 
parameter. We note that if the latter is zero, one has 
to take into account the next-to-leading order bound- 
ary spin relaxation terms, which should determine the 
observable spin density. The solution of the differential 
equation IL'21 12: il) is straightforward and we find for the 
transparent boundary (e.g., ballistic contacts) 



Px{x) 



1 



Re 



k 2 + l 



Re 



ke 



k 2 + 1 



(24) 



and 



Pz(x) 



Im (e zkx ) 



2Re 



(-*- 

\k 2 + l 



(25) 



2£ 

Px{x) = ^-rRc 



and 



Pz(x) 



Refc 



-Im 



ke lkx 1 



fc 2 - 1 



Refc 



fc 2 + l 

fc 2 " - T 



(26) 



(27) 



where k = ^J(l + iV7)/2 » 0.978 + 0.676i is an eigen- 
value of the diffusion equation l|22l \'2:>li . Since the cor- 
responding differential operator is non-Hermitian, the 
eigenvalues are complex. This leads to oscillations of 
the spin density near the boundary with the period 
l osc Ri 6.422V Dt s . We should point out that there defi- 
nitely exists other types of small-lengthscale oscillations 
(with the period tt/pf), which are superimposed onto the 
regular behavior (|24l - 127(1 . However, these Friedel oscil- 
lations are beyond the reach of the diffusion approxima- 
tion. From Eqs. (|24l - [2*7 | . we see that the spin density 
rapidly decays as we move away from the boundary. The 
corresponding lengthscale is l s ~ 1A79\ZDt s . 



V. SUMMARY 

In this paper, we have developed a general method of 
deriving boundary conditions for spin diffusion equation. 
The method involves a gradient expansion of the bound- 
ary Green's function, which takes into account the be- 
havior of the electron wave-functions on ballistic length- 
scales near the edges. Using the proposed method, we 
have found a general form of boundary conditions in the 
Rashba model for two types of edges: a transparent in- 
terface and a hard wall. 

We have found an exact solution of the diffusion equa- 
tion satisfying these two types of boundary conditions. 
The spin accumulation Q24I - I27|l is shown to oscillate 
and decay away from the boundary with the correspond- 
ing lengthscales being determined by the eigenvalues of 
the bulk diffusion equation. However, we argue that even 
small changes in the boundary potential would lead to a 
different solution of the diffusion equation: The eigen- 
values, which determine the density profile, are robust 
and universal, but the overall amplitude and phase shifts 
are boundary-specific. Therefore, the spin Hall effect is 
generally a non-universal phenomenon, which depends on 
the structure of the sample edges. 
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